Boundary conditions for Density Gradient corrections in 3D Monte Carlo simulations

Monte Carlo remains an effective simulations methodology for the study of MOSFET devices well into the decananometre regime as it captures non-equilibrium and quasi-ballistic transport. The inclusion of quantum corrections further extends the usefulness of this technique without adding significant computational cost. In this paper we examine the impact of boundary conditions at the Ohmic contacts when Density Gradient based quantum corrections are implemented in a 3D Monte Carlo simulator. We show that Neumann boundary conditions lead to more stable and physically correct simulation results compared to the traditional use of Dirichlet boundary conditions.